Quantum spin hydrodynamics and a new spin-current mode in ferromagnetic metals 
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We derive the quantum spin hydrodynamic equations in a ferromagnetic metal. From these equations 
we show the existence of a new massive spin-current mode. This mode can be observed in neutron 
scattering experiments and we discuss the difficulties in seeing it. At the end we discuss the existence 
of this mode in localized ferromagnets. 
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Propagating collective modes in condensed matter sys- 
tems are often a consequence of spontaneously broken 
symmetry. In ferromagnetic systems, the spin rotational 
symmetry is spontaneously broken leading to the exis- 
tence of spin waves. They were first observed in iron us- 
ing neutron scattering experiments in the nineteen fifties 
Theoretically spin waves were predicted |||| in lat- 
tice models with a ferromagnetic ground state as well 
as in ferromagnetic Fermi liquids 0. In paramagnetic 
metals, spin waves propagate in an external magnetic 
field || . Not all propagating modes are a consequence of 
spontaneously broken symmetry and an example is zero 
sound in liquid helium Q . 

In this paper we will demonstrate the existence of a 
new collective spin mode in ferromagnetic metals. We 
start from a quasi-classical kinetic equation, the Landau 
kinetic equation, for the momentum distribution func- 
tion in a ferromagnetic metal. From this we can de- 
rive the dynamic equations for the evolution of the spin 
density and spin current. In the low temperature limit 
the spin dynamics is collisionless and dominated by spin 
precession, and we will call this the quantum spin hy- 
drodynamic (QSH) regime. At higher temperature this 
dynamics crosses over to the classical spin hydrodynam- 
ics (CSH) of Halperin and Hohenberg @]. In the QSH 
regime we show that a new massive mode exists along 
with the usual magnon. While gaped, it exists outside 
of the particle-hole, Stoner, continuum and therefore it 
can propagate. From a macroscopic point of view, this 
mode is induced by the collective oscillations of the spin 
current and it should be observable in neutron scattering 
experiments. 

In the long wave length limit one can derive the ex- 
istence of the spin waves from the conservation law j7]] 
for the magnetization in the CSH limit. However, the 
spin current is not a conserved quantity and therefore 
one cannot derive propagating modes related to the fluc- 
tuations of the spin current from the conservation laws. 
Here the Landau kinetic equation first used by Abrikosov 
and Dzyaloshinksii Q for weak ferromagnetic metals is 
used to derive the full QSH theory. 

We consider a three dimensional, weak itinerant fcrro- 



magnet below its Curie temperature. To avoid compli- 
cations we ignore the effects related to the lattice. The 
phenomena which we will describe is independent of the 
lattice, although lattice effects can add more features to 
the ones that we study. We assume that the itinerant fer- 
romagnet is well described by ferromagnetic Fermi liquid 
theory . The kinetic equations for the magnetization 
and the spin current can be derived from the equation for 
the Fourier transform of the Wigner distribution function 
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with ap + ^ 2 ,/3 the single fermion operator for a particle 

with spin (3. Here Dq — ^nyi an< ^ e p^{^i *) i s the effec- 
tive quasiparticle Hamiltonian. Also Iq is the quantum 
collision integral which contains contributions from two 
body collisions and from the rate of change due to the 
spin precession M. 

Following Ref. || we expand all single body operators 
in the basis set spanned by the Pauli matrices 
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at every point (r, t) of space and time with ap the mag- 
netization and 
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the effective magnetic field which is a sum of an external 
magnetic field Ti. and an internal magnetic field gener- 
ated by the quasiparticle interactions /2-, . We expand 
the quasiparticle interactions in spherical harmonics and 
separate the spin symmetric (/ ; s ) and spin antisymmetric 
(/") parts. The result is 
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Here we treat the quasiparticle interaction as spin rota- 
tion invariant. This is possible in the case of small mag- 
netizations and external magnetic fields as explained in 
pot - Here N(0) is the average density of states, averaged 
over the two Fermi surfaces. We must note that the con- 
dition that the ground state of the metal is ferromagnetic 
means that iV(0)/ o Q = F a < -1. 

For an inhomogeneous system the spin quantization 
axis is position dependent and therefore the spin density 



o(r,t) = 2 J Dpop{f,t) 
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is position dependent. From the kinetic equation one can 
obtain the continuity equation for the spin density which 
reduces to the usual spin conservation law 
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where the spin current is 
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and a summation over repeated indices is used. The pre- 
cession term in the r.h.s. of Eq.(10) follows from the non- 
commutativity of the particle density operator and the 
quasiparticle Hamiltonian. Assuming that the magneti- 
zation is small and that we are close to equilibrium we 
can linearize the spin current which leads to 



js,i{r,t) = 2 Dpvp.op 1 
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because in the Fermi liquid only momenta in the vicinity 
of the Fermi surface contributes to the magnetization. 
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is the equilibrium distribution func- 



tion, and F" is a Fermi liquid parameter. 

Starting from the kinetic equation for the spin density 
distribution function 
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then multiplying by the quasiparticle velocity and inte- 
grating over the quasiparticle momentum the linearized 
equation for the spin current takes the form 



where the spin stress tensor is 
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Next we apply a small, transverse (to the external 
field and therefore to the magnetization) magnetic field 
SH = 5TL x e x + STL y e y which induces a change in the mag- 



netization Sap 
vectors and Sap 



3ap x c x 



Sapyiy, where e$ are the unit 
o^. To make the expansion in 



spherical harmonics more transparent we define 
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where Soi — Sop x ± iSapy, o^ is the magnitude of the 
magnetization in the absence of the perturbing field, and 
m is the equilibrium magnetization. We expand in spher- 
ical harmonics the quantities 



(18) 



and substitute the corresponding expressions in Eq.(13). 
Here k is the unit vector in the direction of the spin 
quantization axis. Now Eq.(13) involves terms grouped 
according to the value of the integer number I. Next 
we neglect the terms with I > 1. Although the term 
I = 2 is very important as we discuss later, the derivation 
which we present below is clearer and contains all the 
steps needed for the derivation of the dispersion including 
higher spherical harmonics. Then using the expansion in 
spherical harmonics, Eq.(18), and in the presence of an 
external magnetic field, Eq.(14) reduces to 
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Here we have used a relaxation time approximation for 
the collision contribution to the spin current. 

Before we examine the full spin dynamics coming from 
Eqs.(lO) and (19) it is useful to contrast this with the 
behavior of a paramagnetic metal. In a paramagnetic 
Fermi liquid the first two terms on the r.h.s. of Eq.(19) 
are proportional to col — 277i, where Ti. is the external 
magnetic field. At high temperature the CSH regime is 
realized and at low temperature the system crosses over 
to the precession dominated regime. The cross over scale 
is set by lu^td ks 1. This behavior is similar to that in 
paramagnetic metals In the ferromagnetic metal 

the cross over scale is nearly independent of the external 
magnetic field. To see this we first note that at high 
temperature, Eq.(19) in the uniform case reduces to 
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which means that jg , exponentially decays and its ex- 
plicit time dependence can be ignored. In this limit, 
which is called the classical spin hydrodynamic limit with 

H = the equations reduce to those used by Halperin 
and Hohenberg B and it is easy to derive the dispersion, 
lo r~j q 2 for the spin waves B. The quadratic depen- 
dence for the Goldstone mode follows from the fact that 
the order parameter m is a conserved charge of the spin 
rotation algebra. As the system crosses over to low tem- 
perature it enters into the QSH regime in which the spin 
current changes in time due to the precession around the 
effective field with cu^ = 2m(/f/3 — f§ ). For zero exter- 
nal field the cross over from QSH to CSH occurs when 
uj \ t d ~ 1- The range of temperatures where oj^tjj >> 1 
is what we have referred to as the QSH regime. In this 
QSH regime we can not ignore the -§[jff,i term in Eq.(19). 
This as we will see is the source of the new spin-current 
mode. 

There is another way to see the difference between 
spin waves in a ferromagnetic and a paramagnetic metal. 
Consider a ferromagnetic metal in an external magnetic 
field. The magnetization of the metal will align with the 
external field. If we apply a small pulse in a direction not 
collinear with the external magnetic field, the magnetiza- 
tion will start to precess around the external field. Next 
we reduce to zero the external magnetic field. In a para- 
magnetic metal the polarization will also go to zero and 
the spin wave will become damped. However, in a ferro- 
magnetic metal the magnetization will continue to pre- 
cess about the same axis, although no external magnetic 
field is present to define that axis. Now one can imagine 
looking at the system from a frame of reference rotating 
with the magnetization. In this frame of reference there 
will be no magnetization precession, but there will be a 
precession of the spin current, because the spin current 
is not a conserved quantity and as seen from Eqn.(19) it 
precesses with the frequency u>i . 

It is instructive to look first at the modes of the system 
in the simpler case when Fj 1 = 0, for I > 1 and the ex- 



ternal field is TL = 0. In that case the equations simplify 
to 

dSa(r,t) d . ,„„, 
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where c 2 = |1 + Ftf | (1 + -3-)-^- is the spin wave veloc- 
ity, with Vf the Fermi velocity. Here we have introduced 
S3 = c?(r, t) — <t°, with 3° = mk , where m is the magneti- 
zation. If the system is in the QSH regime and 83 -C 3° 
one can write for the spin current the following hydrody- 
namic equation 
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Differentiating this equation with respect to Xi and using 
the equation for the spin density we obtain 
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Below we will set h = 1. Applying again a transverse 
magnetic field as above Eq.(16) we can write 
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Assuming that the solution is of the form 

5a+ = m+e 1 ^-^ (26) 
with q = qe x we obtain the dispersion equation 



u 2 + c 2 sq 2 + 2 m (f*- f -±)u = 
The two modes have the dispersion 
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where the I = 0(+) solution corresponds to the spin 
waves and the I = 1(— ) solution corresponds to a new 
collective spin mode. For small q we have the usual Gold- 
stone mode 



^o(g) 

and the new mode has the dispersion 
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Under close examination one finds that to the q 2 term 
there is a contribution from the I — 2 spherical harmon- 
ics in Eq.(18). Taking that into account gives the full 
expression for this mode 
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The above equation is valid for qvF << Tn\ff |. In the 
limit /" — > the mode merges with the Stoner continuum 
and is Landau damped. For /" ^ the mode propagates 
with a dispersion that depends on the sign of /". 

It is useful to look at the contribution of this mode 
to the /-sum rule. It is known that the spin waves do 
not satisfy this sum rule. However the spin waves plus 
the new mode exhaust the /-sum rule and the spectral 
weight is shifted from the Stoner excitations to this new 
mode. To see this we first introduce the dynamic struc- 
ture function 
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where x(g, u>) is the dynamic spin susceptibility obtained 
from the kinetic equation by noting that x — jm- Here 
0(uj) is the step function. From the expression for the 
dynamic structure function, calculated from the kinetic 
equation Eq.(13), we can derive the f-sum rule for the 
transverse dynamic spin response function, (q,ui). 
This is given by 



1 F? nn 2 

^S+-(q>) = (l + ^-)^L, (33) 



where n is the total density and m* is the effective mass. 
This has the same form as the one obtained for the spin 
dynamic structure function in the paramagnetic phase 

This sum rule is not exhausted by the Goldstone mode 
alone. If we now include the new mode at zero tempera- 
ture we have 



7+- 



(q, u) = a+-S(uj - u (q)) + P+~ S(oj - Wl (g)) (34) 



and this exhausts the /-sum rule, Eq.(33). This is sig- 
nificant, because the spectral weight of the Stoner exci- 
tations has been transfered to the two modes, leaving 
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This might be a plausible explanation for the difficulty 
of the direct observation of the Stoner excitations in a 
neutron scattering experiments Jl2[ in Fe and Ni since 
the oscillator strength of these single excitations has been 
reduced. 

The observation of the new mode described in this pa- 
per is difficult because it has a small spectral weight at 
small q. This follows from the q 2 dependence of 
term in Eq. (34) . Another obstacle in observing the mode 



is that it originates from the oscillations of the spin cur- 
rent, which is not a conserved quantity and therefore it 
is easily damped. 

In conclusion, in this paper we described a new collec- 
tive mode in weak ferromagnetic metals in the quantum 
spin hydrodynamic regime. While our calculations are 
for small moment itinerant ferromagnets this mode will 
exist in other ferromagnetic metals in which the moment 
is not small. In that case the calculations are significantly 
more involved. This new mode should exist as well in lo- 
cal moment ferromagnets. Mathematically this can be 
seen if one relaxes the condition used by Halperin and 
Hohenberg pj] that the partial time derivative of the spin 
current is zero. This as well as a calculation of this mode 
at finite temperature we leave for a future publication. 
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